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ABSTRACT 

Broad-band images in the Ca II H line, from the BFI instrument on the 
Hinode spacecraft, show emission from spicules emerging from and visible right 
down to the observed limb. Surprisingly, little absorption of spicule light is 
seen along their lengths. We present formal solutions to the transfer equation 
for given (ad-hoc) source functions, including a stratified chromosphere from 
which spicules emanate. The model parameters are broadly compatible with 
earlier studies of spicules. The visibility of Ca II spicules down to the limb in 
Hinode data seems to require that spicule emission be Doppler shifted relative 
to the stratified atmosphere, either by supersonic turbulent or organized spicular 
motion. The non-spicule component of the chromosphere is almost invisible in the 
broad band BFI data, but we predict that it will be clearly visible in high spectral 
resolution data. Broad band Ca II H limb images give the false impression that 
the chromosphere is dominated by spicules. Our analysis serves as a reminder 
that the absence of a signature can be as significant as its presence. 

Subject headings: Sun: chromosphere 
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Introduction 



The Hinode spacecraft is a stable pl atform from which unique high resolution, seeing- 
free images of the Sun can be acquired flKosugi et a/.l 120071). The BFI instrument, fed by 
the Solar Optical Telescope (SOT) on Hinode (ITsuneta et a/.ll2008l ). can observe a 3 A wide 
spectral bandpass centered at the H line of Ca II. Over this bandpass, the line forms in both 
the photosphere in the wings, and chromosphere in the core. Movies of such Ca II images 
have revealed a remarkably dynamic, spicule-dominated limb. The observed spicules have 
smaller diameters, higher app arent velociti e s and smaller lifetimes ( Ide Pontieu et a/.l 120071 ) 
than was previously thought ( iBeckerd Il968l . Il972l ) . 



Figure [T] shows a typical snapshot from a series of Ca II BFI images acquired on 7 
November 2007, in the northern polar coronal hole. We selected a coronal hole because 
spicules there are longer than elsewhere, thereby providing a broad background of spicule 
emission against which a stratified atmosphere might easily be identified. A smooth radial 
gradient has been divided out of the data to enhance the emission high above the limb. The 
zero point of the height scale z = (along a; = 0) corresponds to the sta ndard fo r matio n 
height of the 5000 A continuum, when observed vertically, as derived by lBj0lsethl ( l2008l ). 
Henceforth we will refer to heights on this standard scale. Bj0lseth found that the Ca II 
limb lies 0.45 ±0.034 Mm above the blue limb. Since the continuum at the limb forms about 
0.375 Mm higher than at disk center. 



le.g. 



Athayi 119761 ). the Ca II limb forms near heights 



of 0.825 Mm. 



Curiously, such limb images show little or no signature of a bulk, stratified chromo- 
spher^ which, as we discuss below, should have a thickness between 1 and 2 Mm. The 
BFI instrument's resolution (~ O'.'l) is ample to resolve structure on scales of 1-2 Mm 
(O'.'l = 0.0725 Mm). Yet a striking feature of these images is the continuous emission 
seen along each spicule all the way down to, and sometimes across, the limb. Where then 
is the stratified chromosphere, and why are spicules so obviously dominant that one might 
conclude that the chromosphere itself consists of little more than a collection of spicules? 
In this paper we explore what these observations imply in terms of the structure of the 
chromosphere. 



^By "chromosphere" we refer not only to the traditional definition of Ha emitting plasma seen during 
eclipse flashes, but all the material lying between the quiet photosphere, with densities above ~ 5 x 10^^ 
g cm~^, and the corona with densities below ^ 10"^'^ g cm~^. 
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2. Simple calculations 

The complex dynamic behavior of the Ca II spicules, their unknown origin and other 
difficulties preclude the possibility of meaningful ab-initio or other detailed modeling efforts. 
To address the above questions, a far simpler calculation is appropriate. We model the 
chromosphere as a stratified atmosphere from which spicules emanate. Formal solutions to 
the transfer equation along rays tangential to the solar limb are performed for prescribed 
source functions, densities and atomic parameters applicable to the Ca II H line. 



Our adopted stratified medium is simply the ID atmosphere of iGingerich et al\ ( 1l97ll ). 
As suggested by the Hinode BFI observations, this medium does not emit as much as the 
embedded spicules, so it serves primarily to scatter photons. The exact stratification of the 
bulk chromosphere is not critical, all that is required is that it span the range of densities 
from photosphere t o corona in < 2 Mm . The mean stratifications of hyd r odyna mic or MHD 



models, such as by ICarlsson and SteinI (119951 ): IWedemeyer-Bohm et all (120091 ). are similar 



to the stratification used here. 

The stratified atmosphere was assigned source functions of By{T{r))f{r) where T{r) 
is taken from the atmospheric model, as a function of radial distance from Sun center r. 
By(T) is the Planck function at frequency z/ and temperature T. The function!^ /(r) mimics 
well-known non-LTE effects in which source functions fall below their LTE values. Figure [2] 
shows the source functions used in our calculations. 

The spicules were treated statistically, both in their spatial distributions and thermo- 
dynamic properties. They were randomly distributed along the boundaries of circular su- 
pergranules into 200 small bushes with a common "root", 8 spicules in each bush. 1600 
spicules per supergranule, each with a diameter of 0.1 Mm, leads to a ffiling factor by area 
of 0.015, and a total of 2 x 10^ such spicules on the Sun at any time. (The latter is some 20 
times larger than the value derived by Beckers 1972, based upon data of far lower angular 
resolution than those from the Hinode BFI instrument). These numbers produce synthetic 
spicule images similar to those from Hinode. 

The roots of individual spicules are set at tq = 1.25 Mm above the continuum pho- 
tosphere, with a base density p = po = 1-4 x 10^^^ g cm~^. Below the tq there are no 
spicules in our calculations- the atmosphere is 100% "stratified" . The spicules are modeled 
simply as longer extensions of the atmosphere from this base. The spicule densities are 
p(r) = po 6"'-'""''°^/^= with i.s = 3.5 Mm for r > tq. The scale height £ was chosen to match 



^Here f{r) — 1 (r < ri) and e ''O/i 'i'^ where ri = tq + 4.4£, £ = 0.125 Mm being a typical pressure 
scale height in the photosphere. 
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intensity scale heights of ~ 3.5 Mm found for polar coronal holes by (lBj0lsethll2OO8l ). in order 
to compare calculations with Figure [1] The calculated intensities depend only weakly on the 
gas densities, because our source functions are fixed and the spicules have optical depths of 
order 10^ in the line cores. The assumed spicule properties will be revisited in the Section HJ 

Spicule orientations were randomized relative to the local vertical in azimuth, and their 
source functions specified along each one's length but randomly varied between spicules. 
The source functions are not individually known, being determined locally by colhsional 



excitation and scattering of radiation from the bright underlying photosphere. iMakital (12003 



fig. 12) shows source functions below 4 Mm with black body temperatures near 4300K. Here, 
each spicule's base source function was chosen from a randomly distributed sample about a 
mean of By{T = 4300K) with an arbitrary distribution width one tenth of this. Along each 
spicule the source function drops with height along with the density. Fig. |5] shows the mean 
value as a function of height. For the line opacity, all calculations use a calcium logarithmic 
abundance of 6.3 (where H=12), all calcium is assumed to be in Ca II (a good approximation 
below the transition region) an absorption oscillator strength of 0.33, a microturbulence of 
10 km s~^ (except where specified below) and a radiat ively damped Voigt profile. Standard 
continuum opacity was added to the line opacity from lAUeru (119731 ) . 



3. Results 



Figure |3] shows intensity profiles of the H line as a function of wavelength and height 
in a "standard" model. The computed intensities are similar to observations both above 
and below the limb. The chosen Ca II line parameters prod uce results not incompatible 
with typical profiles seen on the disk ( Linskv and Avrett 197oh . and with Ca II observations 



obtained both during and outside eclipses ( BeckersI 



1968 



1972tlMakita||2003|). 



Figure H] shows emergent intensities at various monochromatic wavelengths and inte- 
grated over the BFI filter bandpass. In the "standard" calculation (left panels), the spicule 
line profiles are assumed to be the same as in the stratified atmosphere. The latter leads to 
absorption at wavelengths within 0.1 A of line center and below heights of ~ 2 Mm. The 
path lengths and line opacity of the stratified material are sufficient to produce absorption, 
unless a spicule happens to lie physically closer to the observer. The limb in the simulated 
BFI data (see the panel labeled "H inode BFI Ca 11") is close to the 0.875 Mm value derived 
observationally by lBi0lsethl ( 120081 ). In the lowest panels, the same radial function was ap- 
plied to the simulated data as the observations shown in Figure [H to enhance the visibility 
of spicule emission over the limb. In the broad BFI bandpass, a significant and observable 
fraction of spicule emission is absorbed by the stratified atmosphere below 2-3 Mm. This 
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behavior is inconsistent with the appearance of Hinode data. 



Real spicules are dynamic, as seen both throu gh linewidths and physical motions (IBeckers 



1968. 


1972; 


Makita 


2003: 



Makital l2003t Ide Pontieu et al\ |2007| ) . Therefore we made two further calcula- 



tions: one using broad spicule emission line profiles, and another using spicule-aligned flows. 
Both calculations shift the spicule emission outside of the absorption profiles of the stratified 
atmosphere when the Doppler speeds exceed tq, where £ is the chromo spheric micro- 



turbulence (< Cs, the sound speed, ~ 7 kms ^, e.g., IVernazza et a/.lll98ll ). and Tq > 1 is 



the line center optical depth tangential to the limb. For values of tq varying between 10 and 
10^° the required shifts are a few times c^. 

The right hand panels of Figure H] present calculations including a spicule line broadening 
microturbulent parameter drawn from a distribution with a mean of 30 kms~^ and a stan- 
dard deviation of 10 kms~^. This supersonic microt urbulence is c ompatible with linewidths 
measured from eclipse data below about 4 Mm (e.g. iMakital 120031 ). (In the dynamic calcula- 
tion, not shown, spicule-aligned outflows were drawn from a distribution with a mean of 120 
kms~^ and a standard deviation of 40 kms~^.) In this calculation, the spicules can be seen 
down to and crossing the limb, as observed, and the dark absorption band resulting from the 
stratified atmospheric absorption is less pronounced. Figure |5] shows intensities in the BFI 
bandpass averaged along the direction tangential to the limb, normalized to limb values, as 
a function of height, from the three calculations and from observations. The filter- integrated 
emission from broad spicule line profiles is larger than from the standard calculation, because 
the computed individual spicules are optically thick across their axes, at least for heights 
below 4 Mm. 

The differences in the broad dips in intensity between heights of 1 and 2.5 Mm shows 
our main result- dynamical calculations are needed to avoid such a large dip in BFI Ca II 
intensities across the limb. 



4. Discussion 



It seems that the absence of the stratified chromosphere in the images obtained with 
the Hinode BFI Ca II filter may be explained, at least in part, simply by large Doppler 
shifts resulting from spicule dynamics. The computations (Figure S]) resemble observations 
(Figure [T]) when spicule emission is Doppler shifted out of the dark core of the H line in the 
stratified chromosphere. 
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4.1. Prom photosphere to corona 



The solar atmosphere does not end at the visible photosphere- there must exist material 
as the upwardly stratified extension of the photosphere. This material must, on av erage, be 



highly stratified because quiet Sun coronal pressures are close to 0.1 dyncm ^ (e.g. iMariska 



19921 ). yet photospheric pressures are orders of magnitude higher. The only question of 
interest here is if this material is expected to be able to scatter the Ca II resonance lines. 
Even if the chromosphere were in hydrostatic and radiative equilibrium, and hence maximally 
stratified, almost all of the calcium would be in the Ca II ground state, and the stratified layer 
would span ~ 1 Mm before coronal conditions were reached. In semi-empirical ID models the 
transition from photosphere to corona spans 1 .5 Mm (measured fro m temperature minimum 



to corona. 



Gingerich et a/.lll97ll : lAthay 



19761 : IVernazza et al\\l98l\ \. 



This transition must, on average, be stratified close to hydrostatic equilibrium, because 
motions observed in spectral lines formed in the photosphere and chromosphere are, statis- 
tically speaking, sub-sonic. Indeed one has to look hard to find the on -disk counterparts of 



the highly supersonic type II spicules (IMcIntosh and De Pontieul 120101 ). for example. More 



directly observable signatures of the stratified chr omospheric m edium are found, for exam 



pie, in the "fiash spectrum" seen during eclipses ( iMakital l2003l . and references therein), or 
in the upward extension of photospheric wave motions seen on the solar disk. While the 
observationally-defined "chromospheric extent" inferred by fiash spectra exceeds hydrostatic 
values, it is also compatible with a hydrostatic stratification in the first 1-1.5 Mm. The 



Gibson 


(1973 


); 


Athav 


(1976 


); 


Judge 



(I2006h . 



4.2. Validity of our results 



The ad-hoc parameters in our calculations clearly limit their usefulness. But our essential 
result- the need to Doppler shift spicule material out of the absorbing stratified chromosphere 
in order to reproduce qualitatively Hinode Ca II data- is relatively insensitive to such details. 
The result simply requires spicules to originate close to the base of the chromosphere, and 
have different source functions and/or opacities from the stratified atmosphere. Given these 
conditions, and spicule lengths which exceed the thickness of the stratified atmosphere, our 
result appears rob ust. Our parti c ular choice of parameter s were taken from observed prop- 
erties discussed by iMakital (120031 ) ; Ide Pontieu et al\ (120071 ) . For a given density, the opacity 
in the Ca II H line follows from atomic data, ion abundance, and thermal and non-thermal 
motions. Average densities of the stratified medium are, as we argued above, approximately 
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in hydrostatic equilibrium. However, our spicule densities and their height dependence were 
chosen simply to produce computations qualitatively similar to the particular coronal hole 
data shown in Figure [H 



Ph ysical considerations suggest that po cannot greatly exceed 10 g cm ^. 



de Pontieu et al 



(120071 ): iMcIntosh and De Pontieul (120101 ) find these spicules to be highly supersonic, ~ 100 
km s~^. Such high speeds require magnetic forces in plasma where the sound speed is cer- 
tainly < 10 km s~^. While the mechanism driving spicules is not known, the characteristic 
Alfven speed va must exceed 100 km s~^. Using an upper limit of 1 kG for field strengths 
in the low chromospheric network (1 kG is characteristic of network photospheric fields), 
Va ~ 100 km s~^, we find p < 10~^ g cm~^. But this is an unrealistically large estimate, 
since chromospheric magnetic fields are weaker due to geometric expansion of network fields, 
and not all of the local magnetic energy is free to be converted to kinetic energy. Strong 
network magnetic fields tend also to be largely unipolar, thus only tangential components 
associated with magnetic shear or with weaker neighboring opposite polarity fields contain 
the free energy. Adopting field strengths nearer to 0.1 kG, as an order of magnitude estimate, 
the observed spicule speeds require p ~ 10~^^ g cm~^, as used above. It is difficult to see 
how this estimate can be significantly larger. 



Spicules in c oronal holes are longer , as seen i n gro und based data (IBeckersI 119721 ) and 
in Hinode data (Ide Pontieu et al\ 120071 : lBj0lsethl l2008l ). Our scale height of 3.5 Mm for 
coronal hole densities and source functions is twice the value derived for the numbers of 



spicul es observed as a function of height for the Sun in general by iBeckerj ( 119721 ). lBj0lseth 
(120081 ) shows in her Fig. 4.10 that the Hinode Ca II data of equatorial regions have scale 
heights closer to 2 Mm. The relationship between the spicules observed by the Hinode BFI 
instrument and earlier work has not yet been clarified. We simply note that our calculations 
are not unrealistic parameterizations of the conditions needed to describe radiative transfer 
in the chromosphere of a coronal hole. 

Our computations are not in q ualitative di s agree ment with the 1.5 Mm wide dip in Ha 
line center intensities discovered by iLoughheadl ( 11969! ) . The cores of Ha and n eutral heliuna 



lines routinely show a dip in intensity surro unded by a shell of emission (e.g. IWhitelll963 
Loughheadl Il969l : jPope and Schoolman! Il975l . and much later work). However, dips seen in 
visible lines of hydrogen and helium may result more from the well-known lack of opacity 
in the low to mid chromosphere, and extra opacity due to fibrils which appear to over-arch 
the stratified chromosphere. Resolving the issue would require detailed calculations of Ha 
and He lines with models taking into account the fibril structure and excitation mechanisms 
populating these excited atomic levels, which are not currently feasible. 



Lastly, there remains an interesting discrepancy between the Hinode BFI data and fiash 
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spectra, in that the 2 Mm Ca II sca le height is twice the median value derived from flash 
spectra of Ca II lines (IMakital l2003l . his table 2, see also our Fig. [5]). Importantly we also 
note that our calculations never remove the off-limb dip entirely, even in the dynamic and 
broad line calculations, yet at least some of the Hinode BFI images appear to show no hint 
of a dip. 



4.3. Implications 

Both observations and simple physical arguments require that spicules be a consequence 
of some plasma or magnetohydrodynamic processes occurring within the chromosphere. 
Spicules cannot arise fully fledged from the photosphere for several reasons, not least because 
photospheric and spicular densities and gas pressures differ by many orders of magnitude. 

Our result suggests that the Hinode BFI Ca II images can be used as speedometers 
in the sense that spicules, when visible down to the limb, must have components that are 
Doppler shifted supersonically, say by more that say 20 km s~^. While the limb chromosphere 
appears in the Hinode BFI Ca II data to be made entirely of spicules, this broad bandpass 
appears to be almost blind to much of the stratifled, inter-spicule chromosphere. The Hinode 
BFI Ca II fllter probably misses populations of the short "type-I" spicules also associated 
with the magnetic network, whose line widths and Doppler shifts are insufficient to avoid 
the absorption by the intervening material. 

In fact, these Hinode data miss the bulk of the mass of the chromosphere, including 
the internetwork. Standard chromospheri c models give on ave rage 0.03 g cm~^ as the total 



chromospheric surface mass density (e.g. IVernazza et a/.lll98ll ). The mass density per unit 



area of the spicules, averaged over the surface, is ~ pif where p, £ and / are their typical 
mass density, length and surface filling factor. Using p ~ 1.4 x 10"^^ g cm~^, i ~ 3.5 x 10^ cm, 
/ = 0.015, we find an average mass density of only 7x 10~^ g cm~^. The spicules observed by 
the Hinode Ca II BFI instrument comprise less than 0.3% of the entire chromospheric mass. 
The energy fiux density needed to support the net work chromosphere against radiation losses 



is estimated to be a few times 10^ erg cm~^ s~^ ( lAnderson and Athayl 119891 ). The enthalpy 
fiux density of individual spicules with speeds of 100 kms~^ is j9f ~ 3 x 10^ erg cm^^ s^^, 
which is thus comparable. Perhaps then these Hinode spicules are intrinsically related to 
the chromospheric heating that is observed? 



- 9 - 



5. Conclusions 

Hinode BFI Ca II images obtained at the solar limb are consistent with the presence of 
the stratified chromosphere when spicular emission is Doppler shifted relative to the stratified 
material. This can be achieved most naturally using broad and/or Doppler shifted spicule 
line profiles of magnitudes compatible with observed motions. The picture presented here 
can be tested directly using very stable spectra at the solar limb, to see for example if the 
behavior modeled in Figure H] is qualitatively correct. This is a very challenging observation 
to make from the ground, but should be possible under conditions of outstandingly good 
seeing and with modern adaptive optics systems. 

The calculations reinforce a commonly known problem regarding broad band spectral 
imagers: one must be very careful taking care of physical effects such as Doppler motions 
which are not spectrally resolved by the instrument. BFI Ca II limb observations are largely 
blind to the bulk of the chromosphere itself. This fact is a sobering reminder that the absence 
of a signature can be as significant as its presence. 
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Fig. 1. — A snapshot from a time scries of Ca II images obtained with the 3 A wide filter 
of the BFI instrument using the SOT on the Hinode spacecraft. The height scale has been 
carefully determined, relative to the vertical continuum (5000 A) optical depth unity, by 
Bj0lseth. 
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Fig. 2. — Prescribed source functions used in the calculations. The spicules were treated 
statistically. 
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Fig. 3. — Line profiles computed from a radial slice of the "standard" transfer calculations, 
simply to demonstrate that the computed spectra are not dissimilar to those observed. The 
self-reversed core in the photosphere (heights near zero) change to emission in the spicular 
material at greater heights. 
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Fig. 4. — Intensities computed at several monochromatic wavelengths and in the Hinode BFI 
passband are shown as a function of position along the limb tangential diretion and radial 
height. The "standard" spicule conditions were applied (left panel), and broad spicular 
emission lines were computed (right panel). 
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Fig. 5. — Average intensities, normalized to the same point within the solar limb, are plotted 
for observations and for thr ee mode l calcu lations. Also shown is the height dependence of 
the Ca II emission found by iMakital (120031 ) from eclipse flash spectra, the absolute value of 
which is arbitrary. 



